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ABSTRACT 

We present an angular power spectrum analysis of the 2 Micron All Sky Survey 
(2MASS) full release extended source catalogue. The main sample used includes 
518 576 galaxies below an extinction-corrected magnitude of Kg =13.5 and limited 
to \b\ > 20°. The power spectrum results provide an estimate of the galaxy density 
fluctuations at extremely large scales, r< 1000 Mpc. We compare this with mock 
predictions constructed from the ACDM Hubble Volume mock catalogue. We find that 
over the range 1^ I ^100 the 2MASS C/ is steeper than that for the Hubble Volume 
model. However, in the linear regime (/ ^30) there is good agreement between the two. 
We investigate in detail the effects of possible sources of systematic error. Converting 
linear power spectrum predictions for the form of the three-dimensional matter power 
spectrum, P{k), and assuming a flat CDM cosmology, a primordial ns=l spectrum 
and negligible neutrino mass, we perform fits to the galaxy angular power spectrum at 
large linear scales {I ^30, corresponding to r> 50 h^^ Mpc). We obtain constraints on 
the galaxy power spectrum shape of Pe// = 0.14 ±0.02, in good agreement with previ- 
ous estimates inferred at smaller scales. We also constrain the galaxy power spectrum 
normalisation to (crgfe/f)^ = 1.36 ± 0.10; in combination with previous constraints on 
(Tg we infer a i^s-band bias of 6x=l-39±0.12. We are also able to provide weak con- 
straints on flmh and flb/^m- These results are based on the usual assumption that the 
errors derived from the Hubble Volume mocks are applicable to all other models. If we 
instead assume that the error is proportional to the Ci amplitude then the constraints 
weaken; for example it becomes more difficult to reject cosmologies with lower Fe//. 

Key words: cosmological parameters - cosmology: observations - large-scale structure 
of the Universe - infrared: galaxies 



1 INTRODUCTION 

The nature of galaxy fluctuations at extremely large scales 
(r< 1000 Mpc) is poorly constrained. Over the last 
decade, large galaxy surveys have constrained the form of 
the galaxy density field to a few hundred Megaparsecs. How- 
ever, the agreement with the concordance model at these 
scales can only be weakly inferred. Indeed, recent evidence 
has suggested that there may be excess power over the ex- 
pected ACDM form to the three-dimens ional power spec- 
trum of matter, P(k), at large scales jFrith et alj l2005l 
l2004l. l2003l: iBussweU etaP l2004l) arising from large inho- 
mogeneities in the local galaxy distribution. 

In recent y ears, large r edshift su r veys o f both galax- 
ies jCole et al.l |2005; Pcr cival et all I2OO1I : IZehavi et all 
|2002^ and QSOs (Outram et alJ 120031) "have determined 
P{k) at relatively small scales. Using the 2dF Galaxy 
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Redshift Survey (2dFGRS), ICole et al.l (l2005l) have con- 
strained the form of galaxy density fluctuations to scales 
of r «300/i~^ Mpc and the associated cosmological param- 
eters to r2„/!.=0.168±0.016 and r2b/f2™=0.185±0.046 (as- 
suming h—0.72). However, determining the power spectrum 
through such redshift surveys suffers from large statistical 
uncertainty at large scales due to the relatively few objects 
available, as well as uncertainties arising from cosmic vari- 
ance due to the relatively small volumes surveyed. 

Using imaging surveys as opposed to redshift surveys 
provides a greater number of objects over larger solid angles. 
With angular power spectrum analysis of such surveys it is 
therefore possible to constrain the form of galaxy fluctua- 
tions to extremely large scales. However, since the clustering 
signal from a particular scale in real space is smeared over 
a range of angular scales, cosmological constraints through 
comparisons with linear theory predictions at smaller scales 
cannot be mad e; the departure from linearity at scales of 
r<40/i"^Mpc iPercival et al.ll200j) affects the clustering 
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Figure 1. A selection of /Cs-band magnitude estimates from the 2MASS full release compared with lLovedavl i2000l) Jf-band photometry. 
In each case the lower panels display the residual. The x = y slope is indicated by a solid line, while the mean offset is indicated in the 
lower panel by a dashed line. This offset (in magnitudes), the best fit slope determined from least squares fits and the rms scatter are 
indicated for each magnitude estimate. We estimate the magnitudes directly from the (a) ii's-band extrapolated and (c) ifs-band fiducial 
elliptical Kron magnitudes, and also from the (b) J-band extrapolated and (d) J-band circular Kron magnitudes colour-corrected to the 
Ks-ha,nd using the J and ii's-band fiducial elliptical Kron magnitudes. 



signal in the angular power spectrum over a wide range of 
scales. Nevertheless at large scales, where this effect is in- 
significant, angular power spectrum analysis represents one 
of the most effective probes of local large-scale structure. 

Previously, the galaxy angular power spectrum has been 
determined for the Sloan Digital Sky Survey Early Data 
Release, the Edinburgh-Durham So uthern Galaxy Cata- 
logue, and a sample of IRAS galaxies ijTe ^rnark et al ] |2002t 
iHuterer. Knox fc Nicholl 1200 it IScharf et al.lll99l respec- 
tively), which along with the recent analyses of redshift 
surveys has constrained the form of galaxy ffuctuations to 
r « 300/i"^Mpc . 

The 2 Micron All Sky Survey (2MASS) has recently 
been completed and provides Ks, H and J-band photom- 
etry for ~ 1.6xl0^ extended sources over t he entire sky to 
Ks>13.5 Jjarrettil200i: IJarrett et alj|2000l) : at the time of 
writing this dataset represents the largest all sky galaxy sur- 
vey. The 2MASS data therefore represents a uniquely power- 
ful probe of the local galaxy density field at large scales; ap- 
plying a galactic latitude cut of |6| >20° in order to remove 
regions of high extinction and stellar contamination yields a 
sample containing 518 576 Ks < 13.5 galaxies, probing a vol- 
ume approximately 5 times larger than the final 2dF Galaxy 
Redshift Survey (2dFGRS) volume. A further advantage of 
2MASS over previous datasets is that the photometry is 
extremely accurate with high completeness; the photomet- 
ric zero-point calibration is accurate to 2-3 per cent; galaxy 



identification is «99 per cent reliable and the galaxy cata- 
logue is >90 per cent complete for |fe| >20° (Jarrctt 2004). 

In this paper, we use data from the 2MASS final release 
extended source catalogue to determine the ifa-band galaxy 
angular power spectrum with the aim of determining the 
form of the clustering of galaxies at extremely large scales, 
and constraining the shape and normalisation of the power 
spectrum. In section 2, we describe the 2MASS dataset and 
the magnitude estimator used. The method of analysis is 
outlined and the 2MASS angular power spectrum is deter- 
mined and compared to mock power spectra in section 3. In 
section 4, we investigate various sources of systematic error. 
We determine constraints for various cosmological parame- 
ters in section 5. The conclusions follow in section 6. 



2 DATA 

2.1 The 2MASS Extended Source Catalogue 

The 2 Micron AU Sky Survey (2MASS) final release ex- 
tended source catalogue provides Ks, H and J-band pho- 
tometry for over 1.6x10^ extended sources over the entire 
sky with high completeness to A"s = 13.5 Ijarrett 2004V 

Previously, in order to estimate the total Ks-hanA mag- 
nitudes from the 2MASS second incremental release data, 
ICole et all (120011 ) used the deeper J-band Kron magnitudes, 
colour-corrected to the Ks-haxiA via the J and Ks default 
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aperture magnitudes. The accuracy of this magnitude es- 
timator was determine d through a co mparison with the 
Tf-band photometry of iLovedavl l|200 rf): the Loveday pho- 
tometry had better signal-to-noise and resolution than the 
2MASS scans and so enabled more accurate 2MASS magni- 
tudes to be determined. 

The final release data uses revised magnitud e estimates 
and t he default aperture magnitudes used in ICole et all 
l|200lfl have been abandoned (Jarrett - priv. comm.). In 
Fig. 0we show a selection of 2MASS Ks-hanA magnitude 
estimate s wit h the revised 2MASS photometry compared 
with the lLove dav (2000) photometry used previo usly. In the 
place of the default aperture magnitudes used in lCole et alJ 
i200ll) . we use fiducial elliptical Kron magnitudes in pan- 
els (b) and (d) to colour-correct the J-band magnitudes to 
the T^'a-band. Of the many different magnitude estimates 
examined, the most accurate in terms of the scale error be- 
tween the Loveday and 2MASS photometry and the zero- 
point offset uses the J-band extrapolated magnitude colour- 
correcte d to the JCg-band as d escribed above. Using the dust 
maps of ISchleeel et alj (Il998f) , the main galaxy sample uses 
extinction-corrected /('s-band magnitudes calculated in this 
way. 

In order to verify the usefulness of the magnitude esti- 
mator used in this work as an estimate of the total Ks-hanA 
magnitude, we perform an internal check via a comparison 
with the magnitude estimates used in the 2MASS-selected 
6dF Galaxy Survey (6dFGS). The 6dFGS if.-band magni- 
tudes are determined using a surface brightness correction to 
the i^s-ban d 20 mag, arcsec" ^ isophotal elliptical aperture 
magnitude i Jones et al.ll20o3l . We find excellent agreement 
with a slope of 1.022, an offset of 0.018 magnitudes and a 
spread of cr=0.048 magnitudes for |6| >20° galaxies matched 
below Ks = 13.5. 

The 2MASS dataset removes or flags sources identified 
as artefacts such a s diffraction spikes and meteor streaks 
J Jarrett et alj|2000l) : we use the 2MASS cc.flag to remove 
such objects. We also employ a colour cut (J — Ks <0.7 
and J — Ks >1.4) below Ks=12 in order to remove a small 
number of objects identifie d asnon-extragalactic extended 
sources iMaller et alJ l2003l l2005^ . In this work, our main 
sample includes 518 576 Ks <13.5 galaxies above a galactic 
latitude of |&| =20°. For reference, the surface density is 19.1 
deg~^. We also use a shallower sample limited at Ks — 12.5 
and |6| >20° which includes 124 264 galaxies and for which 
the surface density is 4.58 deg~^. 



2.2 The ACDM Hubble Volume Simulation 

The Hubble Volume catalogues represent the largest volume 
N-body simulations of the Universe to date. The ACDM 
simulation follows the evolution of 10** dark matter particles 
from z ~50 over a volume of 3000^ h~^Mpc^ to a resolution 
of ~3h~^ Mpc. The associated cos mological parameter s are 
nm=0.3, D.b=OM, h=0.7, (78=0.9 l|Jenkins et al.lll99 ?^. 

In this work, we construct mock 2MASS catalogues 
from the z — ACDM Hubble Volume simulation dark 
matter particles. We divide the total volume into 27 vir- 
tually independent spherical volumes of r = 500 Mpc. 
These are subjected to the 2MASS selection function: 
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jBaugh fc Efstathioul Il993l: iMaller et all l2005ll where z is 
determined from the 2MASS-2dFGRS matched sample de- 
scribed in Frith et al. (2005); for reference 2=0.074 for 
Ks <13.5 and 2=0.050 for Ks <12.5. Equation [T] is nor- 
malised to match the total number of observed 2MASS 
galaxies for |b| >20°. Due to the volume of the 27 mock 
2MASS catalogues, the selection function is artificially trun- 
cated for the Ks <13.5 mocks at z ~0.156. However, this has 
a negligible effect on the work in this paper; at this redshift, 
~95 per cent of the galaxies are sampled for Ks <13.5. 



3 THE 2MASS ANGULAR POWER 
SPECTRUM 

3.1 Estimating the Power Spectrum 



Following the usu al m eth od (e.g. [ Peebl ej ^73 
Hauser fc Peebled 119731 iPeebles fc HauseJ Il974 



^^^^^^^^^^^^^^^^^^^^^^^^^^^ ^^^^^^^^^^^^^^^^^^^^^^^^^^^^ ^^^^^^^^^ 

ScIw^^jT^9^ . the angular power is estimated through 



a spherical harmonic expansion of the surface density of 
galaxies. The coefficients of this expansion are determined 
over the observed solid angle blobs'- 



Yr{e,(t))dQ. 



(2) 



where M=Ngai/^obs is the observed number of galaxies per 
steradian. The angular power is then determined: 
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m 

where, 

JT= / vneA 



(3) 



(4) 



(5) 



The angular power is then normalised, subtracting the ex- 
pected shot noise contribution: 

such that C;=0 corresponds to a random distribution. 

3.2 Fitting to the Power Spectrum 

In order to compare the angular power spectrum with cos- 
mological predictions, we determine an expected form for the 
angular power spectrum for various cosmological parameters 
using the relation between the three and two-dimensional 
power spectra: 



I m|2 2 

\ai \ = - 

TV 



r'^^{r)ji(kr)dr] k^P(k)dk+N' 



(6) 



JScharf et al.l Il992l: iTeemark et all l2002l: 

IlL^ere^^^nox^ NicholT 1200 j) . which we normalise as 
before. Here, $(r) is the 2MASS selection function, and j; is 
a spherical Bessel function. The 2MASS selection function 
is determined using equation 1. 

We use the tran sfer function fitting formulae of 
lEisenstein fc Hul il998l) to obtain a linear theory prediction 
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Figure 2. The |b| >20° 2MASS galaxy angular power spectrum for (a) 518 576 Ks < 13.5 and (b) 124 264 Ks < 12.5 galaxies. The 
crosses indicate the 2MASS datapoints with the shaded region and solid line indicating the Icr spread and mean power spectrum of 
the 27 mock unbiased 2MASS catalogues constructed from the ACDM Hubble Volume mock catalogue as described in section 2.2. In 
each case, unbiased linear theory models corresponding to the Hubble Volume mock catalogue input parameters of f2m=0.3, f2i,=0.04, 
h=0.7 and as=0.9 are indicated by the dashed lines. In the lower panels we show the fractional deviation of the 2MASS power spectrum 
from this model applying the best fit power spectrum normalisation, asb^=1.3(i, (determined in section 5 for the Kg <13.5 sample) to 
the linear prediction, with errors taken from the mock 2MASS la spread. In addition we indicate the approximate mean distance scale 
probed by the data for each i-mode on the top x-axis. 



for the dark matter power spectrum, P{k), with input pa- 
rameters for the matter, vacuum, baryon and neutrino den- 
sities (57„i, 57a, and Q.^,), h (such that Hq — 100/i) and 
matter power spectrum normalisation (as). We also employ 
a hnear biasing scheme such that Pgai{k)=b'^ Pmatter{k) to 
provide a linear prediction for the galaxy P{k). This is then 
transformed to a galaxy angular power spectrum prediction 
using the spherical Bessel function transform in equation |S| 

3.3 Results 

The angular power spectrum for 518 576 Ks < 13.5, |6| >20° 
2MASS galaxies is presented in Fig. determined through 
a spherical harmonic expansion of the galaxy number den- 
sity as described in section 3.1. In order to determine the 
expected scatter due to cosmic variance we determine the 
angular power spectrum for the 27 unbiased mock 2MASS 
catalogues constructed from the ACDM Hubble Volume sim- 
ulation described in section 2.2; the mean and la spread are 
indicated by the solid line and shaded region. On the top x- 
axis we also indicate the approximate distance scale probed 
by the angular power spectrum at the mean depth of the 
sample determined from the 27 mock 2MASS catalogues. 
At the very smallest Z-modes, the Ks <13.5 power spectrum 
probes scales of > 500 Mpc. 



We have also calculated the linear prediction corre- 
sponding to the ACDM Hubble Volume input parameters 
{^^=0.3, nA=0.7, ^6=0. 04, h^O.7, a8=0.9 and 0^=0) 
through a spherical Bessel function transform of the three- 
dimensional power spectrum as described in section 3.2; this 
is indicated for a bias of 1.0 by the dashed line. The linear 
model and the mean mock 2MASS power spectrum are in 
good agreement at large scales. At smaller angular scales 
{I >30) the effects of non-linear clustering become signifi- 
cant. 

In order to verify whether the form and scatter of the 
mock power spectra, which we later use to estimate the er- 
ror on the observed angular power spectrum, is consistent 
with the data, we perform a fit between the two. We 
marginalise over the normalisation of the mean mock angu- 
lar power spectrum and use the binning as shown in order 
to reduce the covariance to insignificant levels. We assume 
that the spread in the mock power spectra is independent 
of normalisation, i.e. we apply the same spread determined 
for the unbiased mock power spectra to the observed an- 
gular power spectrum. In this particular case, this is likely 
to provide an optimistic view of the observed errors since 
we are not shot noise limited. In this scenario, the errors 
are likely to be independent of the power spectrum ampli- 
tude; on the other hand, if the observed power spectrum 
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Figure 3. The \b\ >20°, Ks <13.5 extinction-corrected 2MASS galaxy angular power spectra for the four magnitude estimators shown 
in FiE^. [ Tl using (a) the raw magnitude estimate and (b) a zero-point correction to account for the offset determined with respect to the 
lLoveda^l2000l) photometry. The la spread and mean 2MASS mock power spectrum are shown as in Fig.|21 The lower panels indicate the 
effect of each magnitude estimator on the resulting power spectrum compared to the colour-corrected J-band extrapolated magnitude 
estimator (with the zero-point correction) used in Fig.l^Jand also indicated here by the black crosses in panel (b). In the upper panels 
we have displaced the kmext and kmfo datapoints for clarity. 



is cosmic varia nce limited the errors s cale with model nor- 
malisation (see iFeldman. Kaiser fc Peacock (1994 ) for fur- 
ther discussion on this point). We investigate the impact of 
this assumption on the associated cosmological constraints 
in section 5. First, we perform a fit over the full angu- 
lar range 1< I sJlOO between the Ks <13.5 2MASS galaxy 
angular power spectrum and the mean mock 2MASS power 
spectrum; we find that / d.o.f .—'i.Q. Limiting the angu- 
lar range to scales which are not significantly affected by 
non- linear clustering (Z ^30), the form of the mock power 
spectra are in better agreement with the observed 2MASS 
galaxy angular power spectrum, with /d-o.f. =2.0. 



The form of the 2MASS angular power spectrum is 
therefore in good agreement with the ACDM prediction in 
the linear regime, although it is clear from Fig.|2^ that there 
is some difference in slope at small scales. Assuming the va- 
lidity of the prediction, this is due either to scale-dependent 
bias in the non-linear regime or resolution effects in the Hub- 
ble Volume simulation. Consistency with the ACDM predic- 
tion in the linear regime, of interest in this work, is confirmed 
through a comparison (in the lower panel) with the linear 
prediction for the ACDM Hubble Volume simulation input 
parameters applying a scale-independent bias to match the 
normalisation of the observed power spectrum at large scales 
(see section 5). 



4 SYSTEMATIC ERRORS 
4.1 Magnitude Limits 

Before turning to the cosmological constraints inferred from 
the 2MASS galaxy angular power spectrum it is important 
to verify that the results are robust and not significantly af- 
fected by potential sources of systematic error. While the 
2MASS catalogue is >98 per cent reliable for |fe| >20°, 
Ks <13.5 galaxies jjarrett et al.l200oll and 99 per cent com- 
plete for l&l >30°, 12. 0< Ks <13.7 galaxies jMaller et alJ 
I2OO3), we wish to verify that the angular power spectrum 
is robust to changes in the magnitude limit, and is not ad- 
versely affected by variable incompleteness or reliability at 
faint magnitudes or scale errors in the photometry. 

Fig.|21shows the 2MASS galaxy angular power spectrum 
as a function of imposed magnitude limit. The shape and 
normalisation of the power spectrum, with respect to both 
the linear model and the mean mock 2MASS power spec- 
trum, are remarkably robust to changes in the magnitude 
limit. The departure of the linear model from the observed 
power spectrum occurs at larger angular scales with the shal- 
lower magnitude limit due to the reduced mean depth of 
the sample. For this reason also, the mock 2MASS power 
spectrum is more significantly distorted at the very small- 
est scales by resolution affects resulting in a slightly steeper 
slope at l> 70. 
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Figure 4. The |b| >20°, Ks <13.5 2MASS galaxy angular power 
spectra including no extinction correc tion, and as previous ly an 
extinction correction derived from the lSchlegel et al] Jl99^ dust 
maps. The mock 2MASS mean angular power spectrum and la 
spread are shown as before. In the lower panel we indicate the 
effect of this correction on the power spectrum through a com- 
parison with the corrected sample (indicated by the crosses in the 
upper panel and as shown in Fig|^). 



4.2 Magnitude Estimator 

Throughout this paper, we estimate the Ks-hand magni- 
tudes using the J-band extrapolated magnitudes colour- 
corrected using the Ks and J-band fiducial elliptical Kron 
magnitudes, as this results in a smaller zero-point offset and 
scale error when compared to the more accurate K-hand 
photometry of Lovedav (20i]3)- We wish to investigate the 
effect on the power spectrum by the choice of magnitude 
estimator; in Fig. |2K and b we compare the power spectra 
for the four magnitude estimators presented in Fig. Q with 
and without respectively the associated correction to the 
Loveday ( 200fl) zero-point. 

The power spectrum is robust to changes in the mag- 
nitude estimate and zero-point at the < 10 per cent level. 
This is due to the fact that the change in the depth of the 
survey due to differences in the magnitude limit and scale 
error effects are insignificant. 

4.3 Extinction 

While the level of extinction in the Ks-hand is low and 
the 2MASS magnitudes have been corrected using the 
FSchlegcl ct al... (^1998) dust maps, it is useful to examine the 
potential level of systematic error introduced by extinction. 
Fig. 2] shows the 2MASS galaxy angular power spectrum 
with and without correction for extinction. In this extreme 
case, the effect of removing the dust correction to the mag- 



Figure 5. The |6| >0°, Kb <13.5 mean power spectrum and la 
spread determined from the 27 mock 2MASS catalogues (solid 
line and shaded region). As in Fig. |^ the dashed line indicates 
the expected linear trend for the ACDM Hubble Volume mock 
input parameters of f2m=0.3, 05=0.04, h=0.7 and O"8=0.9. As 
a consistency check, we also show the ACDM Hubble Volume 
non-linear power spectrum (large dot-dashed line) calculated via 
the numerically-determined ACDM Hubble Volume P{k) (Carl- 
ton Baugh - priv. comm.) transformed to the angular power spec- 
trum as described in section 3.2. In the lower panel we compare 
this prediction with the \b\ >0° (dot-dashed line) and |6| >20° 
(solid line) mean mock 2MASS power spectra. 



nitude estimate is at the < 10 per cent level at large scales 
and < I per cent above I ~10. 



4.4 The Window Function 

Throughout this paper a |ti| >20° galactic latitude cut is ap- 
plied in order to avoid the high levels of extinction and stellar 
contamination in the zone of avoidance. We wish to deter- 
mine the level of any systematic effect on the spread of the 
Hubble Volume mock power spectra (and therefore our in- 
terpretation of the statistical uncertainty) introduced by the 
window function. In Fig. |^the 27 mock 2MASS power spec- 
tra and corresponding linear theory model for the ACDM 
Hubble Volume input parameters are shown with no galactic 
latitude cut. Neither the shape nor the spread of the power 
spectra are significantly altered. The effect of the window 
function on the angular power spectrum is < 5 per cent at 
all scales. 

In order to check the consistency of our results we pro- 
vide a further verification of the mock 2MASS power spec- 
trum results through a comparison with the transform of the 
numerically-determined ACDM Hubble Volume simulation 
P(fe) (Carlton Baugh - priv. comm.). There is good agree- 
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Figure 6. Filled contours representing the Itr, 2a and 3(T confi- 
dence regions for the galaxy power spectrum shape and normali- 
sation determined from fits to the 2MASS |6| >20° <13.5 
galaxy angular power spectrum in the range I ^30. The cross indi- 
cates the best fit parameters of rj,j'j^=0.14 and ((786/^)^=1.36. We 
also show the Icr confidence region for the 2MASS result as above 
where we use errors which scale with the model power spectrum 
normalisation (dashed line). 



ment with both the j 
power spectra. 



>0° and Ifel >20° mean mock 2MASS 



Figure 7. Contours of decreasing likelihood in the Qm/i - r2j,/Qm 
plane for the best-fitting angular power spectrum in the range 
I S530. The filled contours indicate the Icr, 2ct and Scr confidence 
regions for the 2MASS |6| >20° Ks <13.5 galaxy angular power 
spectrum, determined from simple fits, marginalising over the 
normalisation and h. The solid contours indicate the Itr, 2(t, 3(t 
and 4(T confi dence regions determ ined from the 2dFGRS 100k 
release P{k) <Percival et"ai]|200ll) . The cross marks the best fit 
model to the 2MASS data of O.mh=0A2 and ^5/0^=0.39. As 
in Fig. |S]we also show the Icr confidence region for the 2MASS 
result as above where we use errors which scale with the model 
power spectrum normalisation (dashed line). 



5 COSMOLOGICAL CONSTRAINTS 

Using the 2MASS galaxy angular power spectrum we have 
determined the form of the galaxy density field at extremely 
large scales and verified that it is not significantly affected 
by common sources of systematic error. We now wish to 
determine the associated cosmo l ogical constraints. 

Using the lEisenstein fc Hul il998ft transfer function fit- 
ting formulae we have determined linear theory predictions 
for the three-dimensional power spectrum of matter, P{k), 
using input parameters of Q.m , , ^Ib, h and matter power 
spectrum normalisation, era; in the subsequent analysis we 
assume a negligible neutrino mass density, a primordial 
ns = l spectrum and i}A=l-^im- We form galaxy angular 
power spectrum predictions using the spherical Bessel func- 
tion transform described in section 3.2 and a linear biasing 
scheme. 

First, we perform fits to the galaxy power spectrum 
shape and normalisation. Assuming a CDM cosmology, the 
power spectrum can be defined through a parameterisation 
of the shape 

Fe// = iimh expi-Qbil +V2h/nm)) (7) 
iSugivamal Il995f) . and a normalisation, which for galaxy 



power spectra may be parameterised through the galaxy 
bias and erg. We constrain Fe// and (crgfcx)^ using a 
grid of 200x800 models between 0.1< F^// sJO.S and 
0.0^(o-8&k)^ <8.0 respectively. We perform least squares 
fits to the \b\ >20°, Ka <13.5 angular power spectrum as 
shown in Fig. |2K at scales of Z ^30 (binned as shown to re- 
duce the covariance to insignificant levels); beyond / ~30 the 
angular power spectrum begins to be significantly affected 
by non-linear effects. 

We take the spread determined from the 27 mock 
2MASS angular power spectra in order to estimate the 
errors on the 2MASS datapoints, assuming that the un- 
certainty remains the same for a biased as for an unbiased 
distribution (as in section 3.3). In doing this, we assume 
that the ACDM Hubble Volume mock catalogue provides 
an accurate description of the local galaxy distribution 
at large scales and that the associated uncertainty in the 
datapoints is realistic. However, since these errors are valid 
only in an unbiased ACDM cosmology we are required to 
make assumptions as to the nature of the cosmic variance in 
the various other cosmologies scrutinised in these fits. Here 
we assume that the errors are independent of cosmology 
and power spectrum normalisation; the likely impact of this 
assumption is examined below. We find that: 
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2MASS clustering dipole of = 1.37 ± 0.34 jMaller et alJ 
1^003) and the 2MASS angular correlation function analy- 
sis of fex = 1.1 ± 0.1 ijMaller et alJl200 5t). The constraint 
on the bias derived in this work rejects 6_ff = 1 at > 3(t; 
it appears therefore that galaxies selected in the T^s-band 
are clustered more strongly than both the underlying mass 
distribution and galaxies selected in optical waveba nds for 
which & «1 fe.g. lVerde et al.ll2002l: rGaztanagalll994ri . 

We are also able to provide constraints on other cosmo- 
logical parameters. We fit to Q,mh and Q.b/^m since these 
primarily determine the shape of the input P{k) and the 
size of the baryon oscillations. We determine model angular 
power spectra in a 71x51x11 grid between 0.1^ ^mh ^50. 9, 
0.0s; fifc/ri™ ^0.5 and 0.4^ h ^0.9 (the effect of h on 
the angular power spectrum is fairly small and we there- 
fore use a lower resolution) . We perform least squares fits to 
the Ks <13.5, \h\ >20° angular power spectrum at scales of 
I ^30, using errors determined for the 2MASS datapoints as 
before which are independent of power spectrum normalisa- 
tion. 



Figure 8. The angular power spectrum for \b\ >20° Ks <13.5 
2MASS galaxies (as in Fig. is compared to a linear theory 
SCDM prediction using input parameters of r2m=1.0, S7i,=0.04, 
/i=0.50 (dot-dashed line), a ACDM prediction using the Hubble 
Volume input parameters as before (dashed line), and the best fit 
power spectrum shape (for I ^30) of re//=0.14 (dotted line). In 
each case we use the best fit normalisation of (a-gb/f)^=1.36. The 
errorbars indicate the Icr spread determined from the 27 mock 
2MASS power spectra. In the lower panel we show the fractional 
deviation from the best fit T f=0.1A prediction. 



Te// = 0.14 ±0.02 
and 

(crgfeA-)^ = 1.36 ±0.10 

marginalising over the normalisation and power spec- 
trum shape respectively. The associated confidence regions 
are indicated by the filled contours in Fig. |S| 

This val ue of Teff is in exce llent agreement with the 
2dFGRS fit JPercival et alj 1200 j) of F^ff = 0.18 ± .04 
(for ;i=0.7) and the WMAP value llSpergel et al J 120031 of 
Fe// = 0.15 ± 0.01 (for W3=0.99). H o wever , our value is 
slightly higher than the iMaller et alJ ll2005ll result which 
constrains Teff ~ 0.116 ± 0.009 at 95 per cent confidence 
using a measure of the three-dimensional A's-band galaxy 
power spectrum via an inversion of the 2MASS angular cor- 
relation function. 

Our constraint on the A's-band galaxy power spectrum 
normalisation of ((Ts^y)^ = 1.36±0.10 is also slightly higher 
than the .Mailer et al. (2005J result of crgfeif = 1.0±0.1. Using 
the WMAP-2dFGRS best fit matter power spect rum nor- 
malisation of as = 0.84±0.04 l|Bennett et alJl200,'^ ). we con- 
strain the Ks-hand bias to bK = 1.39 ± 0.12, in reasonable 
agreement with previous measurements determined from the 



The filled contours in Fig. |7| show the associated con- 
fidence regions marginalising over the normalisation. We 
are able to provide weak constraints on the cosmology of 
Qmh <0.62 and Slt/f^m <0.46 (at la confidence). These 
constraints are particularly insensitive to the baryon den- 
sity since the acoustic oscillations detected in redshift survey 
analyses are smoothed over a wide range of angular scales. 
However our constraints are in good agreement with the 
previous results at smaller scales from the 2dFGRS P{k) 
ilPercival et"ai]|200ll: ICole et al.ll2005^ . As an example of 
how our results can differentiate between different cosmo- 
logical models we show the 2MASS galaxy angular power 
spectrum compared with ACDM and SCDM predictions in 
Fig.lHl 

We also wish to examine our assumption, used through- 
out this work, that the uncertainty due to cosmic variance 
determined from the 27 ACDM mock 2MASS catalogues is 
independent of the power spectrum normalisation. To do 
this, we instead assume that the errors determined from the 
ACDM mock catalogues simply scale with the model power 
spectrum normalisation as would be the case in the cos- 
mic variance limited scenario, and compare the two cases. 
In Figs. Inland |7| we show the associated la confidence re- 
gions by the dashed lines, marginalising over the power 
spectrum normalisation. We find best fit parameters of 
Fe//=0.125±0.030, {asbKf=1.47toA7, 0.b/Qm <0.52 and 
f2m/i < 0.71. This constraint on the galaxy power spectrum 
normalisation implies a Ks-hsLiid bias of 6x=1.441q'^4 (us- 
ing the WMAP-2dFGRS constraint on as as before). It is 
clear that while the associated confidence regions for each 
parameter are slightly larger the results are in fair agree- 
ment whichever error analysis is used. However, it is clear 
from Fig.|S]that using this alternative assumption about the 
errors it is more difficult to reject combinations of high bias 
and steeper Fe// slopes. For example, Fe//=0.05 would only 
be rejected at 2.5o". More simulations of other cosmologies 
are needed to check whether these errors or the errors used 
elsewhere in this paper are most likely to be correct. 
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6 CONCLUSIONS 

We have used 518 576 Ks < 13.5, |6| > 20° galaxies selected 
from the 2MASS full release extended source catalogue to 
determine the associated angular power spectrum and con- 
strain the form of galaxy fluctuations to Gigaparsec scales. 
We have compared this to a ACDM N-body mock predic- 
tion constructed from the Hubble Volume simulation; it is 
in reasonable agreement although there is a discrepancy in 
the slopes at / >30 in that the 2MASS result is significantly 
steeper than the mock prediction. We compare these to a 
linear theory prediction using the ACDM Hubble Volume 
simulation input parameters; there is good agreement with 
the mock prediction at scales where non-linear effects are 
insigniflcant (I ^30). 

Possible sources of systematic error were investigated. 
We first examined the effect of imposed magnitude limit; 
the 2MASS angular power spectrum slope was robust with 
respect to the 2MASS mock and model predictions. The 
2MASS galaxy angular power spectrum is also robust to 
different magnitude estimators and zero-point corrections 
(imposed to agree with the Lovcdav (2000) photometry) at 
the ~10 per c e nt lev el. We correct for extinction using the 
ISchleeel et al.l Jl998l) dust maps; the effect on the angular 
power spectrum is «10 per cent at l< 10, and ~1 per cent at 
smaller scales. Our results are also robust to window func- 
tion effects; the effect of a |6| > 20° cut is < 5 per cent at 
all scales. 

Finally, we have used linear theory predictions for the 
2MASS galaxy angular power spe ctrum formed from the 
transfer function fitting formulae of lEisenstein fc Hul lll998ll 
to determine constraints on f2,„/i and Q.b/^m assuming a 
flat CDM cosmology, a primordial ns = l spectrum and a 
negligible neutrino mass. Our results arc in agreement with 
the 2dFGRS P{k) constraints (Pcrcival ct al. 2001), and 
we are able to provide weak constraints of Q.mh <0.62 
and Q.h/^m <0.46 (at la confldence). Angular power spec- 
trum analysis is particularly insensitive to the baryon den- 
sity since any associated baryon oscillations are likely to 
be smoothed over a wide range of angular scales. How- 
ever, given the huge volume probed (Ri5 times the final 
2dFGRS volume) the associated constraints on the power 
spectrum shape and normalisation are more significant. 
We also determine constraints for the galaxy power spec- 
trum shape, Fe//, and normalisation, (o-g&if)^. In agree- 
ment with the 2dFGRS and WMAP values, we find that 
Fe// = 0.14 ± 0.02. This is slight l y high er than an alterna- 
tive value found bv iMaller et alJ ^2005^ using the 2MASS 
dataset of Fe// = 0.116 ± 0.009, determined through an in- 
version of the angular correlation function. We also tightly 
constrain the T^s-band galaxy power spectrum normalisa- 
tion to [ashKf = 1.36 ± 0.10. Using the WM AP-2dFGRS 
value of as = 0.84 ± 0.04 jBennett et alj|2003l) . this imphes 
a Tfs-band bias of 6k = 1.39 ± 0.12. 

We also investigated the likely impact on our assump- 
tion that the errors which we use to constrain various cos- 
mological parameters, determined from the unbiased ACDM 
mocks, are independent of cosmology and power spectrum 
normalisation by instead assuming that these errors simply 
scale with the power spectrum normalisation as would be the 
case in the cosmic variance limited scenario. We find that 
while the associated confidence regions are slightly larger 



the results are in fair agreement. However it becomes less 
easy to reject models with lower Fe//; therefore although 
the data appears to prefer a ACDM power spectrum slope, 
it may still not be possible to rule out a ssignificantly steeper 
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